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STRUCTURE AND STABILITY OF ONE-DIMENSIONAL 
DETONATIONS IN ETHYLENE- AIR MIXTURES 


S. Yungster and K. Radhakrishnan 
Ohio Aerospace Institute 
Brook Park, Ohio 44142 


Abstract 

The propagation of one-dimensional detonations in eth- 
ylene-air mixtures is investigated numerically by solving the 
one-dimensional Euler equations with detailed finite-rate 
chemistry. The numerical method is based on a second-order 
spatially accurate total-variation-diminishing scheme and a 
point implicit, first-order-accurate, time marching algorithm. 
The ethylene-air combustion is modeled with a 20-species, 
36-step reaction mechanism. A multi-level, dynamically 
adaptive grid is utilized, in order to resolve the structure of 
the detonation. Parametric studies over an equivalence ratio 
range of 0.5 < cf> < 3 for different initial pressures and 
degrees of detonation overdrive demonstrate that the detona- 
tion is unstable for low degrees of overdrive, but the dynam- 
ics of wave propagation varies with fuel-air equivalence 
ratio. For equivalence ratios less than approximately 1 .2 the 
detonation exhibits a short-period oscillatory mode, charac- 
terized by high-frequency, low-amplitude waves. Richer 
mixtures ( (]) > 1.2) exhibit a low-frequency mode that 
includes large fluctuations in the detonation wave speed; that 
is, a ‘galloping’ propagation mode is established. At high 
degrees of overdrive, stable detonation wave propagation is 
obtained. A modified McVey-Toong short-period wave-inter- 
action theory is in excellent agreement with the numerical 
simulations. 

Introduction 

The basic structure of a detonation wave was determined 
independently by Zeldovich, von Neumann, and Doting 
(e.g., Williams 1 ). The ZND model postulates that detonation 
waves have a steady, one -dimensional structure consisting of 
a lead shock wave followed by a reaction zone that is initi- 
ated by shock heating. The heat release in the chemically 
reactive material sustains the lead shock of the detonation 
wave. The minimum sustainable steady detonation speed 
occurs at the Chapman- Jouguet condition, for which the 
combustion products move at sonic speed relative to the lead 
shock. Detonation waves traveling at speeds higher than the 
Chapman-Jouguet velocity are called overdriven detona- 
tions. 

Subsequent experimental and theoretical studies have 
shown that the ZND structure is unstable. The instabilities 


are responsible for the complex detonation systems observed 
in practice, ranging from cellular three-dimensional struc- 
tures to one -dimensional pulsating and ‘galloping’ detona- 
tions. The objective of this paper is to study computationally 
the nonlinear development of pulsating instabilities in 
one-dimensional detonations. 

Much of the previous computational research on the sta- 
bility of one-dimensional gaseous detonations was based on 
the one-dimensional Euler equations for an ideal gas, with a 
one-step, irreversible chemical reaction governed by Arrhe- 
nius kinetics to model the combustion process. For example, 
Fickett & Woodri Abouseif & Toong , and Bourlioux, Majda 
& Roytburd 4 used this simplified model to study the stability 
of overdriven, piston-supported, one-dimensional detona- 
tions using different computational techniques. They used 
the steady, overdriven ZND detonation wave solution as ini- 
tial conditions. Their results indicated that flow instabilities 
could produce various oscillation phenomena, ranging from 
highly uniform to chaotically irregular pulsating detonations. 
The study of Bourlioux et. a! 4 showed in particular the need 
for sufficient grid resolution around the detonation front to 
avoid the appearance of nonphysical numerical artifacts. In 
addition, for highly overdriven detonations they obtained 
steady solutions with their simplified chemistry model. 

These highly simplified chemistry models are convenient 
for mathematical analysis and algorithm development, but 
are of limited applicability. For example, one-step chemistry 
models cannot model the initiation, branching, and recombi- 
nation steps characteristic of the combustion process at the 
detonation front and therefore cannot reproduce the interac- 
tions between the chemical kinetics and fluid dynamics. This 
coupling determines the dynamics of detonation wave propa- 
gation. 

A few one -dimensional detonation wave studies have 
been conducted with more detailed chemistry models. Fick- 
ett, Jacobson & Schott 5 modeled the chain-branching mech- 
anism with modified Arrhenius kinetics, but considered only 
one chemical reaction. Their model could be tuned to cover 
different values for the length of the induction and reaction 
zones. Their work showed that the structure of the pulsating 
detonations (i.e., the frequency of oscillation) is strongly 
affected by the use of a more realistic chemistry model. 
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Short & Quirk 6 and Sanchez et. af have recently reported 
results using generic three-step global branched-chain chem- 
istry models. Their studies were able to reproduce the essen- 
tial dynamics of chain-branching reaction systems. Their 
results indicated that the ratio of the length of the 
chain-branching induction zone to the length of the recombi- 
nation zone plays a major role in the structure and stability 
of detonation waves. However, the three-step reaction mech- 
anisms used in these studies are not representative of any 
particular fuel-oxidant-diluent system. 

Previous computational studies addressing the stability 
of one-dimensional gaseous detonations with a realistic 
chemistry model have been carried out by Sussman 8 , He 9 
and Yungster and Radhakrishnan 1011 . Sussman used a 
detailed hydrogen-air combustion mechanism to investigate 
the oscillation characteristics of a stoichiometric hydro- 
gen-air mixture for different degrees of detonation overdrive. 
His analysis started from the steady, overdriven ZND solu- 
tion and calculated the transition to the oscillatory propaga- 

o 

tion mode. Sussman concluded that the 
longitudinal-oscillation mode depends on the ratio of the 
heat release time to the induction time. He 9 investigated the 
direct initiation of detonations by an energy source using 
both a one-step reaction model and a detailed combustion 
mechanism for hydrogen-oxygen. Yungster and 
Radhakrishnan 10,11 analyzed the structure and stability of 
one-dimensional hydrogen-air detonations for different mix- 
ture equivalence ratios, pressures and degrees of overdrive, 
using a detailed combustion mechanism. Their results 
showed that a modified McVey-Toong short-period 
wave-interaction theory that takes into account the finite 
heat-release time significantly improves agreement between 
theory and numerical simulations over the entire equivalence 
ratio and pressure ranges 11 . 


drive. The motivation for studying detonations in ethyl- 
ene-air mixtures includes the increasing interest in 
developing pulse detonation engines (PDE) fueled by hydro- 
carbons. Ethylene is a gaseous fuel with one of the shortest 
ignition times among hydrocarbons, making it especially 
attractive for PDE application. In fact, several experimental 
studies on PDEs have utilized ethylene as fuel 1415 . More- 
over, ethylene is one of the most important intermediate spe- 
cies in pyrolysis and oxidation of higher hydrocarbons, 
including JP-10 (Li et a/. 16 ). Finally, one -dimensional pul- 
sating detonations have been observed in ballistic experi- 
ments for hydrogen-air mixture only. It is therefore 
important to determine whether these longitudinal instabili- 
ties are unique to hydrogen, or they are present in other fuels 
as well. 

In this study, we consider a tube closed at one end and 
open at the other and compute the development of the deto- 
nation wave from its (direct) initiation to the final establish- 
ment of an oscillating propagation mode. The degree of 
overdrive is determined by the pressure in the driver gas, as 
described later. This approach is significantly different from 
that used in all previous studies 11 ; in particular it eliminates 
the need for the steady-state ZND solution. 

Numerical Simulations 

The conservation form of the one-dimensional unsteady 
Euler equations for a chemically reacting gas mixture con- 
sisting of n s species can be written as 


dQ + dF 

dt dx 


= W, 


( 1 ) 


where 


On the experimental side. Lehr 13 examined the longitudi- 
nal oscillations in shock-induced combustion flows. His bal- 
listic-range experiments consisted of spherical-nosed 
projectiles fired into premixed combustible mixtures. 
Depending on the experimental conditions, steady or 
unsteady flows were obtained. Two types of unsteady modes 
were observed: one consisting of very uniform, low-ampli- 
tude and high-frequency oscillations, and the other consist- 
ing of less regular, high-amplitude and low-frequency 
oscillations. The oscillations observed in these experiments 
are caused by the same longitudinal instabilities associated 
with the one -dimensional detonations considered here. 

The objective of the present study is to extend our previ- 
ous work 10 ’ 11 on the analysis of the structure and stability of 
one-dimensional detonations by studying ethylene-air sys- 
tems using a detailed combustion mechanism, for different 
mixture equivalence ratios, pressures and degrees of over- 


Q = [pp P 2 < •••’?„ . P«. g ] . 


( 2 ) 


p | u, p 2 w, .... p u, p m“ + p, u(e + p) 


iT 


, (3) 


W = [w v w 2 , .... w n ,0, 0] r . (4) 

L n 

i_ 1 p ; . 

the mixture density, u the velocity, p the pressure, and e the 
total energy per unit volume. The terms {w ( } represent the 
rates of production of species from chemical reactions. 


The numerical method used for solving the governing 
equation set is based on Yee’s second-order spatially accu- 
rate total-variation-diminishing (TVD) scheme 17 , and a 
point implicit, first-order, time-accurate marching algorithm. 
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This method is a subset of a more general class of BDF 
methods considered by Yungster and Radhakrishnan . A 
detailed description of the numerical method is given in ref. 
1 1 . Similar first-order, time-accurate point-implicit methods 
have been used by Sussman , Wilson and Sussman , and 
Matsuo and Fuji to study the longitudinal oscillations 
occurring in shock-induced combustion. 

In order to maintain adequate numerical resolution of the 
detonation wave front, without the need to use hundreds of 
thousands of grid points, a multi-level, dynamically adaptive 
grid is utilized. Figure 1 shows a section of the grid at three 
different times as the detonation wave moves from left to 
right. The grid constantly adapts to keep the detonation front 
within the finest grid level. An arbitrary number of levels can 
be specified. Nine or ten grid levels were used in the present 
study, and 100 points were included in the finest grid level. 

The ethylene-air combustion was modeled using the 
short chemical reaction mechanism of Li et a/. 16 , which con- 
sists of 36 reactions among 20 species (see Table 1), with 
nitrogen being treated as an inert (i.e., non-reacting) species. 
This mechanism was designed to yield reasonably accurate 
results (ignition times differed from those of a more detailed 
mechanism by less than 30%) over the range of post-shock 
temperatures between 1000 K and 2500 K, pressures 
between 0.5 and 100 bar and equivalence ratios between 0.5 
and 3. This combustion model was recently successfully 
used by Povinelli and Yungster 21 to study the effects of dis- 
sociation and subsequent recombination on the performance 
of pulse detonation engines. In particular, the equilibrium 
Chapman-Jouguet detonation speed and von-Neumann state 
calculated with this mechanism agreed well with the results 
produced by the chemical equilibrium code CEA . 

Results and Discussion 

There are two methods by which detonation can be 
formed in general: (1) by direct initiation, wherein a strong 
shock wave is generated in the tube (with a charge of solid 
explosive or by using a high pressure reservoir) and (2) by 
transition from a deflagration. In this paper, we consider 
development of detonations with direct initiation. Figure 2 
shows a schematic of the initiation process. A high-pressure, 
high-temperature driver gas, consisting of nitrogen, was used 
in a small region next to the head-end of the tube. When the 
computation is started, a shock wave travels to the right and 
an expansion wave propagates to the left, towards the 
head-end. The shock wave is strong enough to initiate chem- 
ical reactions in the combustible mixture. Depending on the 
mixture properties, the shock wave and the combustion front 
subsequently merge and form a detonation wave. The value 
of the driver pressure, p c i r - lv , will determine the degree of 
overdrive of the resulting detonation wave,/ defined as 


f = d 2 /d 2 cj , 

where D is the actual detonation propagation speed and D C j 
is the theoretical Chapman-Jouguet detonation speed. (Since 
the flow behind an overdriven detonation is subsonic, this 
solution is susceptible to degradation by rarefaction waves 
from the closed end. With the overdriven detonation wave 
established, these rarefaction waves will eventually overtake 
the detonation front and reduce its speed towards the Chap- 
man-Jouguet condition. All cases considered in this paper 
were computed for time intervals much shorter than those 
required to encounter this situation). 

The first case considers a stoichiometric mixture of ethyl- 
ene-air at p 0 = 0.2 bar and Tq = 298 K and a driver pressure 
ratio, r p = Pdri^Po^ °f 150. A grid refinement study for this 
case is presented in figure 3. Numerical results are presented 
for four successively refined grids having a minimum spac- 
ing, Ax min , indicated in each figure. The plots show the vari- 
ation in the detonation speed with time. The initiation 
process is similar to that of hydrogen-air mixtures, and is 
described in detail in Ref. 10. 

It can be observed that the structure of the detonation 
wave changes with grid spacing. For the two coarsest grids 
(figs. 3a and 3b), the detonation reaches a steady propagation 
speed of 1928 and 1932 m/s, respectively, after the initial 
transient phase. With the finest two grids used in this study 
(figs 3c and 3d), the detonation reaches a high-frequency, 
low-amplitude propagation mode. Table 2 summarizes the 
results of the grid refinement study. The final average deto- 
nation speed is essentially independent of the grid spacing. 
The frequency and amplitude of the oscillations are also 
nearly identical for the two finest grids. All subsequent cal- 
culations were carried out on the finest grid. 

The theoretical Chapman-Jouguet detonation speed for 
this mixture, as computed with CEA , is D CJ = 1799.4 m/s. 
Therefore, the results presented in figs. 5 and 6 are actually 
those of an overdriven detonation, with a degree of overdrive 
/= 1-15- 

The high-frequency mode of oscillation can be explained 
by the McVey-Toong wave interaction mechanism . Two 
fundamental processes form the basis of this mechanism: 
(i) when a new reaction front is created, compression waves 
(reaction shocks) are generated that travel upstream and 
downstream from the new reaction front; (ii) when an old 
reaction front is extinguished, it must be accompanied by the 
generation of upstream and downstream rarefaction waves 
that have a strength comparable to the reaction shocks. The 
interaction between these waves and the shock front pro- 
duces the high-frequency oscillations. 


NASA/CR— 2003-212586 


3 



Equivalence ratio effect 

Parametric studies were carried out over an equivalence 
ratio range of 0.5 < <f) < 3 for ethylene-air mixtures at 
p 0 = 0.2 bar and T 0 = 298 K and a driver pressure ratio of 
150. The results are shown in figs. 4-7, and summarized in 
Table 3. 

Figure 4 shows the variation of detonation speed with 
time for several lean mixtures. The structure of the detona- 
tion for all the lean mixtures considered is very similar to 
that observed for the stoichiometric case. The frequency of 
oscillation, ft), decreases monotonically from 0.36 MHz for 
<() = 1.0 to 0.20 MHz for (f) = 0.5 (see Table 3). Note that 
although the driver pressure ratio is fixed at 150, the degree 
of overdrive increases for leaner mixtures. 

Figure 5, which gives the shock pressure as a function of 
time for the same ethylene-air mixtures as in fig. 4, shows 
that the average shock pressure decreases for leaner mix- 
tures, but the normalized oscillation pressure amplitude 
remains nearly constant. The oscillation amplitude is com- 
puted from the peak-to-trough shock pressure normalized by 
the average shock pressure. 

Figures 6 and 7 show the variations in detonation wave 
speed and shock pressure with time for rich mixtures. Note 
that the maximum oscillation frequency is obtained for 
<() = 1.2, and the maximum normalized oscillation pressure 
amplitude is obtained for <\> = 1.6. The degree of overdrive 
remains nearly constant for all mixtures. For <() > 1.2, the 
detonation wave exhibits the long-period mode of instability, 
characterized by low-frequency, high-amplitude oscillations, 
that is, a “galloping” propagation mode. Both the amplitude 
and frequency of oscillation decrease for richer mixtures. 
The frequency of oscillation of the long-period mode is an 
order of magnitude lower than that obtained for the 
short-period mode. Note also that for the <(>=1.6 and 3.0 
cases high-frequency oscillations appear first, but then disap- 
pear once the long-period mode is fully established. 
High-frequency oscillations are also apparent for the two 
richest mixtures during the portion of the long-period cycle 
when the shock pressure is highest. 

Initial Pressure effect 

The effect of initial pressure on detonation wave struc- 
ture was examined by studying its behavior at a higher initial 
pressure (p () = 0.4 bar, 7 (l = 298 K). The results are shown in 
fig. 8 and summarized in Table 4. For the driver pressure 
ratio considered (r p = 150), the degree of overdrive ranges 
from 1.29 for (]) = 0.6 to 1.12 for the two rich mixtures. 

The structure of the detonation wave at the higher pres- 
sure (fig. 8) is similar to that for the p 0 = 0.2 bar case (figs. 4 


and 6), but the frequencies of oscillation are much higher. 
The effect of equivalence ratio on oscillation frequency is 
more pronounced at the higher pressure. For 0 = 1 .0 the fre- 
quency is 0.74 MHz, and it decreases to 0.475 MHz for 
<() = 0.6. For <(>=1.2 the short-period mode appears first, but it 
subsequently switches to the long-period mode at longer 
times. Further mixture enrichment (<(> = 1.4) produces only 
the long-period mode. 

Driver pressure ratio effect 

The effect of driver pressure ratio, r p , on detonation wave 
structure was investigated for stoichiometric ethylene-air 
mixtures at p 0 = 0.2 bar. Figure 9 shows the detonation shock 
pressure as a function of time for several driver pressure 
ratios. 

For the two lowest driver pressure ratios (r p = 70 and 
100) the detonation wave initially exhibits the short-period 
oscillation mode, but it appears to be transitioning to the 
long-period mode at longer times. For driver pressure ratios 
between 150 and 250 (figs. 9c-9e) the detonation wave 
exhibits only the short-period propagation mode during the 
time frame considered. Table 5 summarizes the oscillation 
frequency and amplitude for all the cases shown in fig. 9. 
The frequency of oscillation increases and the amplitude 
decreases with increasing driver pressure ratio. At driver 
pressure ratios of 275 and above, the high-frequency oscilla- 
tions are observed to decay with time; that is, the detonation 
wave is stabilized at high degrees of overdrive. This behavior 
has also been observed by Bourlioux et al 4 using one-step 
Arrhenius kinetics. 

Comparison with Theory 

As discussed previously, the short-period mode, 
described by the McVey-Toong mechanism, results from the 
successive interactions of reaction shocks and expansion 
waves with the lead shock. The resulting oscillation period, 
T mt , is the sum of the ignition delay time, t t , and the time 
required for the reaction shock to overtake the lead shock. 
Implicit in their mechanism is the assumption of negligible 
heat release time relative to the ignition delay time; that is, 
the (single) reaction shock is created immediately upon igni- 
tion. For the one-dimensional detonation flows considered 
here, T mt is given by Alpert & Toong as follows: 

h 

T = (5) 

mt (1 -M s Y 

where M s is the Mach number immediately behind the lead 
shock. Therefore, the frequency of oscillation, to (= 1 IT mt ), is 
inversely proportional to the ignition delay time. 
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Figure 10 shows a comparison of the computed 
short-period frequency of oscillation with that given by 
equation 5, as a function of equivalence ratio for the 
short-period cases discussed previously. The ignition delay 
time is arbitrarily defined as the time required for the tem- 
perature of the gas to increase by 25 K, and was computed 
using the chemical kinetics code LSENS 24 26 . The frequen- 
cies given by the McVey-Toong short-period wave-interac- 
tion theory are in good agreement with the computed 
frequencies for lean mixtures at low pressures. However, the 
frequencies given by equation 5 are much higher for 
near-stoichiometric and rich mixtures, especially at the 
higher pressure. Thus the largest differences between the 
computed and theoretical frequencies of oscillation occur for 
conditions with the shortest ignition delay time. 

The deviation of the McVey-Toong theory (equation 5) 
from the computed frequency has been attributed to neglect 
of the heat release time, an assumption that is not always 
valid since there are many situations in which the heat 
release time is a significant fraction of the total reaction 
time 11 . Modification of the McVey-Toong theory to account 
for the finite heat-release time was proposed in which a mod- 
ified oscillation period, T mmt , is computed as follows 11 : 

t ■ 

T , = + aA t , (6) 

mm l (i— m 1 qmax’ ’ 

where Af qmax is the interval between ignition and maximum 
rate of heat release and a is an empirical constant. In ref. 11, 
which considered hydrogen-air mixtures only, a was esti- 
mated by least-squares minimization to be equal to 0.6. 

The frequencies of oscillation obtained from equation 6 
for the ethylene-air mixtures considered in the present study, 
using the same value of a (= 0.6) as for the hydrogen-air 
mixtures considered in ref. 11 are also plotted in figure 10. 
Note that accounting for the finite heat-release time has 
resulted in significant improvements in the theoretical fre- 
quencies, when compared with the numerical results. In fact, 
the agreement between the modified theory and the computa- 
tions is excellent over the entire equivalence ratio and pres- 
sure ranges. 

Conclusions 

The one-dimensional Euler equations with detailed finite 
rate chemistry were solved numerically to investigate the 
structure and stability of detonation waves in ethylene-air 
mixtures. Parametric studies over a range of fuel-air equiva- 
lence ratios, initial pressures and degrees of detonation over- 
drive indicate that one-dimensional ethylene-air detonation 
waves exhibit unsteady oscillatory propagation modes for 
low degrees of overdrive. For lean to slightly rich ((f) ~ 1.2) 


mixtures the detonation wave propagates in a high-frequency 
low-amplitude mode, as described by the McVey-Toong 
short-period mechanism. For richer mixtures the detonation 
wave switches to a ‘galloping' long-period mode character- 
ized by low-frequency, high-amplitude oscillations. 

High degrees of overdrive stabilized the detonation 
wave: the high-frequency oscillations decayed with time and 
the detonation wave propagated in a stable mode. Calcula- 
tions with very low degrees of overdrive resulted in complex, 
irregular oscillations, possibly transitioning to the 
long-period mode. 

The frequencies of oscillation obtained with the 
McVey-Toong short-period wave-interaction theory are 
much higher than those obtained from the detailed numerical 
simulations, especially for near-stoichiometric mixtures at 
high pressure. This difference is attributed to the implicit 
assumption in the McVey-Toong theory that the heat release 
period is negligibly small compared with the ignition delay 
time. Modification of their theory to account for the finite 
heat-release time resulted in excellent agreement with the 
numerical simulations for all conditions. Exactly the same 
modification has been shown to be effective for hydrogen-air 
mixtures. 
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Table 1: Short Chemical-Kinetic Mechanism for CoH 4 Ignition and Detonation (Li et a/. 16 ). 


Number Reaction 

A* 

n* 

E* 



Hydrogen-Oxygen Chain 



1 

H + 0 2 

-> 

OH + O 

3.52x10'° 

-0.70 

71.4 

2 

OH + 0 

— y 

H + O 2 

1.15x 10 H 

-0.32 

-0.7 

3 

OH + H 2 

-4 

HjO + H 

1.17x10 s 

1.30 

15.2 

4 

H 2 0 + H 

-4 

OH + H 4 

6.72xl0 9 

1.30 

84.6 

5 

0 + H 2 0 

-4 

2 OH 

7.60x10° 

3.84 

53.5 

6 

2 OH 

-> 

O + H 2 0 

2.45x10-' 

3.97 

-19.0 


Hydropcroxyl Formation and Consumption 



7 b 

H + Oj + M 

-4 

H0 2 + M 

6.76X10 19 

-1.40 

0.0 

8 

H0 2 + H 

— ► 

2 OH 

1.70xl0 14 

0.00 

3.7 

9 

HO? + H 

-¥ 

H 2 + 0 2 

4.28xl0 13 

0.00 

5.9 

10 

HO. + OH 

-4 

H 2 O + 0 2 

2.89xl0 13 

0.00 

-2.1 


Hydrogen Peroxide Formation and Consumption 



11 

2 HO . 

-> 

H 2 O 2 + 02 

3.02 xl0 u 

0.00 

5.8 

12' 

H 2 O 2 + M 

-> 

2 OH 4- M 

1.20xl0 17 

0.00 

190.0 


Direct Recombination 



13 b 

H + OH + M 

-> 

H 2 0 + M 

2.20x 10 22 

-2.00 

0.0 

14 b 

H 2 O + M 


H + OH + M 

2.18xl0 33 

-1.93 

499.0 


Carbon Monoxide Reactions 



15 

CO + OH 

-► 

CO 2 + H 

4.40x10° 

1.50 

-3.1 

16 

CO 2 + H 

— ► 

CO + OH 

4.97x10 s 

1.50 

89.7 


Initiation and Fuel Consumption 



17 

C2H4 + o 2 


C2H3 + HO2 

4.22X10 13 

0.00 

241.0 

18 

C 2 H 4 + OH 


C s H 3 + H 5 0 

2.70x10°° 

2.31 

12.4 

19 

C 2 H 4 4 - 0 


CHj + CHO 

2.25x10°° 

2.08 

0.0 

20 

CjH, + 0 

— y 

CH 2 CHO + H 

1.21 xlO 06 

2.08 

0.0 

21 

C2H4 + H 0 a 

-4 

C 2 H 3 + H 2 O 2 

2.23xl0 12 

0.00 

71.9 

22 

C 2 H 4 + H 

-4 

C 2 H 3 + h 2 

2.25 xlO 07 

2.12 

55.9 

23 d 

C 2 H 4 + H 

-4 

C2H5 ko 

1.17x10'° 

0.03 

-28.0 



koo 

1.52X10 13 

0.03 

13.0 


Vinyl, Methyl and Ethyl Consumption 



24 

C 2 H 3 4- H 

-4 

C 2 H 2 + H 2 

3.00 xlO* 3 

0.00 

0.0 

25 

C 2 H 3 4- O 2 

-4 

CH 2 0 + CHO 

1.70xl0 JS 

-5.31 

27.2 

26 

C 2 H 3 4- O 2 

— > 

CH 2 CHO + O 

7.00 xlO 14 

-0.61 

22.0 

27 

CH 3 + 0 2 

-> 

CH 2 O + OH 

3.30x10" 

0.00 

37.4 

28 

ch 3 + 0 

— > 

CH 2 O + H 

8.43xl0 13 

0.00 

0.0 

29 

C 2 H 5 4- Oo 


C 2 H 4 + H0 2 

2.00x10" 

0.00 

20.9 

30 d 

c 2 H 5 

— ► 

C 2 H 4 + H k 0 

1.00x10'° 

0.00 

126.0 



k<X> 

1.30xl0 13 

0.00 

167.0 


Vinoxy, Ketene, Formaldehyde, Formyl and 

Acetylene Consumption 


31 

CH 2 CH 0 

—¥ 

CH 2 CO + H 

1.05 xlO 37 

-7.19 

186.0 

32 

CHjCO + H 

-> 

CHj + CO 

1.11 xl0° 7 

2.00 

8.4 

33 

CH,0 + OH 

-> 

CHO + HjO 

3.90x10'° 

0.89 

1.7 

34' 

CHO + M 

— > 

CO + H + M 

1.86xl0‘ 7 

-1.00 

71.1 

35 

CHO + 0 2 

-> 

CO + HO 2 

3.00x10" 

0.00 

0.0 

36 

C?H 3 + OH 

—> 

CH 2 CO + H 

1.90 xlO 07 

1.70 

4.2 


•Specific reaction-rate constant k — AT n t~ E,RT \ units mol, cm 3 , s, K and kJ/mol. 
b Chaperon efficiencies: CO, 1.9; CO 2 , 3.8; Ha, 2.5; H 2 O, 12.0; others, 1.0. 
c Chaperon efficiencies: Same as b above except H*0, 16.3. 
d Falloff from Hewson and Williams 

•Chaperon efficiencies: CO, 2.5; COa, 2.5; Ha, 1.9; HaO, 12.0; others, 1.0. 


Table 2: Oscillation frequency and amplitude for various grid spacings (/’= 1.15) 


Figure 

Minimum grid spacing 
(cm) 

Average propagation speed 
(m/s) 

Oscillation frequency 
(MHz) 

Oscillation amplitude 
A D/D 

3a 

3.906xl0" 3 

1928 

— 

— 

3b 

1.953xl0" 3 

1932 

— 

— 

3c 

9.766xl0" 4 

1931 

0.36 

0.091 

3d 

4.883xl0" 4 

1929 

0.36 

0.092 
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Table 3: Oscillation frequency and amplitude for various equivalence ratios (p 0 = 0.2 bar; r p = 150) 


Figure 

Equivalence 

D CJ (m/s) 

Degree of 

Oscillation frequency 

Oscillation amplitude 

ratio ((])) 

overdrive (/) 

(CO, MHz) 

(A p/p) 

4 & 5 (f) 

0.5 

1539.6 

1.36 

0.20 

0.19 

4 & 5 (e) 

0.6 

1624.5 

1.27 

0.24 

0.18 

4 & 5 (d) 

0.7 

1687.0 

1.22 

0.28 

0.18 

4 & 5 (c) 

0.8 

1733.8 

1.19 

0.31 

0.20 

4 & 5 (b) 

0.9 

1770.2 

1.16 

0.345 

0.19 

4 & 5 (a) 

1.0 

1799.4 

1.15 

0.36 

0.20 

6 & 7 (a) 

1.2 

1841.5 

1.13 

0.40 

0.17 

6 & 7 (b) 

1.4 

1864.4 

1.11 

0.068 

0.82 

6 & 7 (c) 

1.6 

1869.9 

1.12 

0.058 

1.17 

6 & 7 (d) 

2.0 

1853.0 

1.13 

0.036 

1.01 

6 & 7 (e) 

3.0 

1747.7 

1.12 

0.021 

0.56 


Table 4: Oscillation frequency and amplitude for various equivalence ratios (p 0 = 0.4 bar; r p = 150) 


Figure 

Equivalence 
ratio ((])) 

D CJ (m/s) 

Degree of 
overdrive (/) 

Oscillation frequency 
(CO, MHz) 

8(a) 

0.6 

1629.4 

1.29 

0.475 

8(b) 

0.8 

1745.5 

1.20 

0.63 

8(c) 

1.0 

1814.5 

1.16 

0.74 

8(d) 

1.2 

1857.0 

1.12 

0.80, 0.13 

8(e) 

1.4 

1877.3 

1.12 

0.12 
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Table 5: Oscillation frequency and amplitude for various driver pressure ratios (p () = 0.2 bar; (]) = 1) 


Figure 

Driver pressure 
ratio (r p ) 

Degree of 
overdrive (J) 

Oscillation frequency 
(CO, MHz) 

Oscillation amplitude 
(A p/p) 

9a 

70 

-1.04 

transitioning 

0.849 

9b 

100 

-1.07 

transitioning 

0.306 

9c 

150 

1.15 

0.36 

0.192 

9d 

200 

1.22 

0.48 

0.160 

9e 

250 

1.29 

0.59 

0.109 

9f 

275 

1.32 

0.64 

vanishing 

9g 

300 

1.35 

stable 

— 

9h 

400 

1.45 

stable 

— 
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Fig. 1. Computational grid at three different times. 
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Fig. 2. Schematic of detonation wave initiation. 



(d) 



Fig. 3. Effect of grid spacing on detonation wave solution. 
C 2 H 4 -Air; p 0 = 0.2 bar; <>=1.0; r p = p drh /p 0 = 150; D CJ = 
1799.4 m/s;/= 1.15. 
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Fig. 10. Short-period frequency of oscillation; comparison 
between computed frequencies (present study) and those 
obtained from the original and modified McVey-Toong theo- 
ries. C 2 H 4 -Air; p 0 = 0.2 bar, r p = 150. 



Equivalence ratio, <]> 


Fig. 11. Short-period frequency of oscillation; comparison 
between computed frequencies (present study) and those 
obtained from the original and modified McVey-Toong theo- 
ries. C 2 H 4 -Air; p 0 = 0.4 bar, r p = 150. 
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